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FINAL REPORT 

ACTIVE DAMPING OF MODAL VIBRATIONS BY FORCE APPORTIONING* 

I- INTRODUCTION 

The objectives of this research have been to develop and to evaluate by 
analysis and numerical simulation a method of active structural damping known 
as force apportioning. This method Is based upon the method used In modal 
vibration testing of Isolating modes by multiple-shaker excitation. For 
either application, active excitation or active damping, one chooses a distri- 
bution of as few forces as possible on the structure which will maximally 
affect selected vibration nodes while minimally exciting all other modes. 

NASA sponsorship of this research has fully supported the production 
of (1) a technical paper presented at a major conference and published in the 
proceedings (and presently under review for journal publication) and (2) a 
Master's degree thesis. NASA sponsorship has also partially supported the 
production of a second Master’s degree thesis, which is indirectly related 
to the principal objectives of this research. Descriptions and bibliographical 
references to these reports are given in Section II. 

During the last few months of the project, additional work was completed 
which has not previoipsly been reported in writing. The results of this work, 
which dealt primarily with imperfection sensitivity of the active damping 
method; are summarized in Section III and are presented in more detail in 
appendices. 

A final oral report was presented by the Principal Investigator at NASA 
Langley Research Center on July 18, 1980. 

'fc 

The NASA Technical Officer for this grant is Mr. Robert Miserentino, NASA 
Langley Research Center. 
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”• SUMMARY OF ACCOMPLISHMENTS, AND REPORT REFERENCES 
II. A The Basic Form of Active Damping by Force Apportioning 

This subject has been described In great detail In a Master's degree 
thesis [1] ar«d In a conference technical paper [2]. Reference 2 Is Included 
as Appendix A of this report. 

The paper has also been submitted for Journal publication. In the 
Journal paper. Section V of Ref. 2 will be simplified considerably. It hap- 
pens that the same force apportioning vector can be derived with a much more 
direct approach not requiring the Intermediate Introduction of a condensed 
model. The new derivation Is Included In Appendix B of this report. 

II.B The Dynamic Stiffness Method for Plane Grillages 

This method, which has been discussed In the Semi-Annual Reports for this 
project, Is the subject of a recent Master's degree thesis [3], NASA sponsor- 
ship supported a part of the theoretical development and computer Implementa- 
tion reported In Ref. 3. A technical paper describing this work will soon be 
submitted for publication In the Journal of Sound and Vibration. 

The method provides exact vibration modes for plane grillage structures, 
rather than the approximate modes which result from finite element discretiza- 
tion. We had originally planned to use these exact modes In numerical simulations 
of active damping, such as those described In Section VI of Ref. 2. However, 
as the project termination date approached, it became clear that the additional 
time and effort required to employ the exact modes, as opposed to finite element 
modes, would have no affect on the ultimate conclusions relative to active 
damping. So the dynamic stiffness method was not employed In the study of 

*Bracketed numbers refer to the References section. 
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active damping. However, the dynamic stiffness method Itself has considerable 
Intrinsic value, so our efforts In developing and reporting It have been 
worthwhile. 

Ill- SUMMARY OF STUDIES NOT PREVIOUSLY REPORTED 

III. A. Accuracy of Numerical Simulations of Active Damping 

Since modal series truncation and Runge-Kutta Integration were used In 
the time history simulattons described In Section VI of Ref. 2, the results 
are not numerically exact, In order to check the accuracy of the numerical 
solutions for several cases, we varied the number of modes retained and/or 
the Integration step size from the values stated In Ref. 2. In all cases 
but one, these variations produced Insignificant changes in the results. The 
one exception was a case In which the integration step size was clearly too 
long relative to the periods of the high-frequency modes retained; an oscil- 
latory Instability resulted which we feel was a numerical artifact and not a 
qualitatively correct solution of the equations of motion. 

III.B. Active Damping of Higher- Frequency Modes 

Previous numerical simulations [1,2] have demonstrated only that the 
lowest flexible modes of a structure can be suppressed effectively, but these 
new results show that force apportioning is applicable also to the damping 
of higher-frequency modes. The results are presented in detail in Appendix C. 

III.C. Some Studies of Imperfection Sensitivity 

Previous reports [1,2] have described active damping by force apportioning 
under idealized conditions. Numerical simulations have shown that the method 
can be very effective provided that all necessary modal parameters of the 
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structure are known exactly and that sensors, actuators, and narrow-band filters 
have ideal operating characteristics. But the effectiveness of the method under 
more realistic conditions Is a matter of more practical Interest, so we have 
conducted a limited study of the effects of non-ideal conditions. Specifically, 
we have analyzed the case In which the controlled structure's modal parameters 
are known only approximately and the case of non-ideal frequency filtering. 

I I I . C . 1 Inexact Modal Parameters for the Controlled Structure 

It has been assumed In previous reports [1,2] that the exact modal para- 
meters were known for the purpose of computing the Ideal force apportioning, 
which Is given In Section V of Ref. 2. These new results demonstrate that the 
effectiveness of the active damping method Is not particularly sensitive to 
somewhat Inexact modal parameters, provided that the approximate mode shapes 
are not drastically different from the actual mode shapes. However, inexact 
modal information leads to coupling between the controlled modes (those modes 
denoted by the index s in Eqs. (10) and (15) of Ref. 2), which produces the 
possibility of long-term system Instability. The effects of inexact modal 
Information are described In more detail In Appendix D. 

III.C.2 Non-Ideal Frequency Filtering 

It has been assumed in previous reports [1,2] that the narrow-band fil- 
ters required by the active damping method have ideal operating characteristics, 
namely. Infinitely narrow pass-bands and Infinitely steep rolloff rates, so 
that only the desired frequency components are passed. The new results ac- 
count for filters having finite pass-bands and finite rolloff rates. Numerical 
simulation results show that the method's effectiveness in suppressing 
Initially existing vibrations in the controlled modes Is quite insensitive to 
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typical filter characteristics. However, non-ideal filters do Introduce the 
possibility of long-term system Instability due to coupling between all modes, 
both controlled and uncontrolled. The effects of non-ideal filters are des- 
cribed In more detail In Appendix E. 

IV. SUMMARY OF COMPUTER PROGRAMS P RE PARED 

The computer programs developed for this research are run on the main 
computers of the VPI&SU Computing Center, an IBM 3032 Processor Complex, and 
an IBM System 370 Model 158. The IBM 3032 runs the 0S/VS2 MVS (Multiple 
Virtual Storage) operating system with JES2 (Job Entry Subsystem 2) as job 
scheduler. The IBM 370 runs under the control of Virtual Machine Fac111ty/370 
(VM/370); under VM, the Conversational Monitor System (CMS) provides inter- 
active tjme-sharing support to users at remote terminals. The Computing 
Center provides Fortran subroutines of the International Mathematics and 
Statistics Library (IMSL). Automated plotters used for this research are 
the Versatec 1200 electrostatic plotter (for quick, Inexpensive plots) and the 
Calcomp 1051 drum plotter (for finished Ink plots on vellum, suitable for 
publication). 

All programs developed for this research are coded In Fortran IV and 
compiled by the Fortran H Extended Compiler. With one exception , all programs 
are run with IBM single precision arithmetic. Two types of programs have been 
written: calculation programs and plotting programs. 

The calculation programs generate mass and stiffness matrices for finite 
element models of plane grillage structures, perform modal elgenso lutlons, 

* 

The exception is EIGENR, the locally developed and undocumented eigensolutlon 
subroutine used to calculate modes of small finite element models. EIGENR 
uses IBM double precision arithmetic. 


and do time history numerical simulations of active damping. All calculation 
programs are run from CMS on the IBM 370, The time history numerical simula- 
tion programs* ACTIV for Coulomb- type damping and VIACDA for viscous-type 
damping, make extensive use of the Interactive mode for entry of data and 
Inspection of Intermediate results by the user. 

Two slightly differing Fortran plotting programs, DES2 and DES3, produce 
Versatec and Calcomp plots, respectively, of time history numerical simulations. 
Each program Is controlled by a different set of job control language Instructions. 
The plotting programs run under MVS. 

Documentation for the programs developed In this research Is fairly com- 
plete but Informal, and It will remain In the possession of the Principal 
Investigator. Appendix B of Reference 1 presents flow diagrams describing 
the operation of programs ACTIV and VIACDA prior to Incorporation of the ad- 
ditional capabilities described In Section III.C. 

The programs are operational at this writing and will be kept operational 
until approximately six months after publication of the journal version of 
Ref. 2. At that time, some of the operating programs will be deleted from 
disk storage, and the source programs will be stored as card decks. The 
decks will remain In the possession of the Principal Investigator for possible 
future use. 

In their present form, the calculation programs are restricted to analysis 
of unrestrained plane grillage structures having relatively few degrees of 
freedom. However, with straightforward modifications, the programs could 
handle any type of structure. The calculation programs constitute a stand- 
alone package, except that ACTIV and VIACDA presently use IMSL subroutines 
for performing the matrix inversion Indicated by Eq. (28) of Ref. 2. However, 
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the matrix to be Inverted Is usually quite small (10 x 10 or smaller), so a 
short user-supplied subroutine could be employed with no loss of capability. 

The programs permit time history plotting of all physical quantities that 
a user Is likely to desire. The user may plot the response In any degree of 
freedom, the force produced by any control thruster, and the energy In any 
mode or combination of modes. Figs. 4-7 of Ref. 2 Include examples of each 
of these types of plots. Although all the time history plots inctuded In 
this report (others than those In Appendix A) represent total energy (I.e., 
the energy In all modes), plots of several other physical quantities were 
generated In each case. Some of these other quantities are very Interesting 
(e.g., the energy In only the uncontrolled modes), but they have been ex- 
cluded In order to limit the length of this report. 

V‘ SOME POSSIBLE REFINEMENTS OF THIS RESEARCH 

The topics discussed In this section are not all of equal significance; 
they are presented below in random order. 

The target mode suppression time t^ is defined by Eqs. (7) and (13) of Ref. 
2 as being the time specified for decay of modal coordinate to a prescribed 
level. Comparison of Figs. 4a and 7a of Ref. 2 shows that this definition 
does not lead to comparable rates of energy decay for viscous and Coulomb 
types of damping with the same value of tg. Hindsight strongly suggests that 
defining t^ on the basis of modal energy decay, rather than modal amplitude 
decay, would be a more rational approach. 

In the development of the mathematically ideal force apportioning vector 
in Section V of Ref. 2 (and the revised version in Appendix B, Section 1 of 
this report), the restriction is imposed that the number of control actuators 
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must equal the number of controlled modes, n. « m. This restriction arises 

fl 

from the necessity for Inverting a partition of the modal matrix, which must 
therefore be square. But In actual practice, It will probably be required to 
control more modes than there are available actuators, m > n.. This practical 
requirement might be satisfied, within the context of active damping by force 
apportioning, by calculation of the pseudo- Inverse of a non-square matrix. 

The mathematical details of the pseudo- Inverse calculation are presented In 
Appendix B, Section 2. However, no numerical study has been attempted. 

One difficult problem which Is encountered In active damping by force 
apportioning Is selection of the control thruster locations so as to minimize 
residual energy. We have addressed this problem In an ad hoc , trial -and-error 
fashion, as Is described concisely In Appendix C of this report and Section 
VI of Ref. 2, and more extensively In Ref. 1. However, It Is probable that 
a more efficient approach could be devised on the basis of relative values 
of the Inner products of modal vectors and force apportioning vectors. This 
approach Is suggested by the discussion In the sixth paragraph of Section VI 
of Ref. 2. 

Section III.C and Appendices D and E of this report discuss the possibility 
of system Instability occurring If the force apportioning and/or the filtering 
do not have the mathematically Ideal forms. We have Investigated this pos- 
sibility simply by Integrating the equations of motion and Inspecting the time 
history plots. It Is obvious that a much more sophisticated stability analysis 
could be conducted. Equations (E.2) In Appendix E are the most general equations 
of motion, accounting for non-ideal apportioning and/or non-ideal filtering, 
and they would be the basis for a quantitative stability analysis. 
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VI. CONCLUDING REMARKS 


The concluding remarM In Section VII of Ref, 2 still hold without re- 
vision. In addition, we may comment now on the new results presented in 
Section III and the appendices of this report, These results are fairly con- 
vincing evidence that the effectiveness of active damping by force apportioning 
is not significantly impaired by realistic conditions, as opposed to the 
idealized conditions considered in Ref. 2. 

The study described in this report has been exclusively theoretical and 
computational. It seems reasonable to conclude on the basis of our results 
that active damping by force apportioning appears on paper to be a practically 
feasible method. Further mt* more definitive evaluation of the method re- 
quires an experimental study. 


REFERENCES 


1. J.-F. M. Barthelemy. "Active Damping of Modal Vibrations by Force Apportioning," 
Master's thesis, Dept, of Aerospace and Ocean Engineering, Virginia Poly- 
technic Institute and State University, March 1980. 

2. W. L. Hallauer Jr. and J.-F. M. Barthelemy, "Active Damping of Modal 
Vibrations by Force Apportioning," AIAA Paper No. 80-0806-CP, Proc. of 
the 21st Structures, Structural Dynamics, and Materials Conference, May 
12-14, 1980; Seattle. Submitted for publication to the t burnal of Guidance 
and Control . 

3. R. y. C. Liu, "The Dynamic Stiffness Method for Computing Exact Vibration 
Modes of Plane Grillage Structures," Master's thesis, Dept, of Aerospace 
and Ocean Engineering, Virginia Polytechnic Institute and State University, 

March 1980. 


10 


AIAA Paper No. 80-0806-CP APPENDIX A 

Proceedings of the 2lst Structures, Structural Dynamics, and Materials Conference; 
sponsored by AIAA, ASHE, ASCE, and AHS; May 12-14, 1980; Seattle, 


ACTIVE DAMPING OF MODAL VI8RATI0NS 
BY FORCE APPORTIONING 

W. L. H«H«uer Jr,* inJ J, ‘F. M. 84rLh«l«iny*» 
D^partmont of Acrospact and Octan Engineering 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 24001 


ABSTRACT 

This paper describes with theory and numeri* 
cal simulation a method of active structural damp- 
ing which requires relatively few discrete con- 
trol thrusters suitably positioned on the struc- 
ture. For each vibration mode which is to be 
damped (target mode), a particular apportioning 
of coherently phased control forces is applied 
which strongly affects that mode while minimally 
exciting all other modes. The force apportion- 
ing used is that which would tune a target mode 
if the structure were being shaken in a modal 
vibration test. In contrast to modal testing, 
however, the forces are varied temporally so as 
to dampen, rather than excite, the target mode or 
modes. 



NOTATION 

{.) 

column matrix, vector 

[ 

transpose of [ ], (.) 

C( )“3.(-)®® 

denotes condensation or partition- 
ing relative to both modes and 
degrees of freedom (d.o.f.) 

[( )^3 

denotes partitioning relative to 
modes only 

(y.Y)*Jv 

Ck] 

n X n stiffness matrix 

[ni] 

n X n inertia matrix 

U1 

n X n modal matvix 

u] • 

. (U«]’’)T 

fs 

force apportioning vector for damp' 
ing of a target mode 

ir 

rth column of [.^3, mode shape of 
rth mode 

Q(t) 

n X 1 action (force and moment) 
vector 

g(t) 

n X 1 d.o.f. vector 

i(t) 

vector of normal or modal coordi- 
nates 


, ll. r => s 
|0, rfs 

[m]f^ generalized mass of the rlh mode 
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Ep(t) energy in the rth mode 

K. feedback constant for damping of 

target mode, Eqs. (8) and (l3) 

m number of target modes 

n total number of d.o.f. in struc- 

tural model 


n. number of modes and d.o.f, in con- 

^ ^^dMSfd model 

"ds " “sV^ ' 

Uj. natural frequency of the rth mode 

Pj target^mode viscous decay criterion, 


r index denoting an arbitrary mode 

s Index denoting a target mode 


t time 


target mode suppression time, Eqs. 
* (7) and (13) 

5e„ * initial value of target mode 

normal coordinate 

c, target mode viscous active damping 

factor 


I. INTRODUCTION 


To suppress the vibrations caused by opera- 
tional or environmental conditions, modern vehicle 
and building structures may require additional 
damping beyond that provided by inherent passive 
mechanisms of energy dissipation. The traditional 
source of artificial damping is a passive device, 
but for Some current applications such devices 
are ineffective or are unsatisfactory for other 
reaso.is, such as excessive weight, For these 
applications, additional artificial damping might 
be provided bv an active system, i.e. a system 
which uses sensing, feedback, and control forces 
or moments to oppose the structural motion. 

Herzberg et al. ^ have observed that the well 
established technoTogy of modal vibration testing 
is adaptable to the control of structural vibra- 
tions. This paper describes a method of active 
damping which is just such an adaptation. 

A possible application of this method is for 
the suppression of vibrations in proposed large 
space structures. These structures will be sub- 
jected to sources of disturbance such as propulsive 
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forces or gyroscopic moments to effect maneuvers i 
but the vacuum of space will provide no atmos- 
pheric damping and it seems likely that internal 
structural damping will generally be low. The ' 
potential significance of active damping for large 
space structures has inspired a great deal of 
current research. No literature review will be 
attempted here, but Ref. 2 includes a number of 
very recent papers and much bibliographic Infor- 
mation on the subject. 

Most provious contemporary studies of active 
structural damping have applied advanced control 
theory, e.g. the papers of Balas3 and Meirovitch 
and Oz^. In contrast, the basic form of active 
damping by force apportioning introduced here in- 
volves no control theory beyond the behavior of 
a one-degree-of-freedom damped oscillator. 

II. GENERAL DISCUSSION OF FORCE APPORTIONING 

The traditional application of force appor- 
tioning is for modal vibration testing with two 
or more sinusoidally varying, coherently phased 
mechanical shakers. The shakers are attached to 
suitable points on the structure being tested, 
and their force amplitudes are apportioned in 
Such a way as to tune a selected vibration mode, 
which we call the target mode. Ideally, all of 
the energy is fed into the target mode at its 
natural frequency, and all other modes, which we 
call the residual modes, are untouched. This 
focusing of energy into a single mode is attrib- 
utable to the specific spatial distribution of 
shakers and their relative force amplitudes. 

In reality, perfect tuning is generally not 
possible with a finite number of shakers. The 
practical objective of modal testing with force 
apportioning is to tune with as few shakers as 
possible while still minimizing excitation of 
residual modes so that target mode parameters 
can be accurately measured. Several methods 
have been developed or proposed for calculating 
the apportioning of shaker force amplitudes to 
produce effective mode tuning, e.g. the methods 
of Lewis and WrisleyS, Ashero, Ibanez^, and 
Morosow and AyreS, 

In modal vibration testing, the properly 
apportioned forces feed energy into the target 
mode and stimulate motion by virtue of being in, 
Phase with the velocity of target mode response. 
Conversely, for active damping of existing but 
unwanted vibration in a target mode, the same 
spatial distribution and apportioning of forces, 
but 180° out of Phase with the velocity, will ex- 
tract energy from that mode. The target mode or 
modes can always be suppressed in this manner, 
but a possible undesirable byproduct is the 
pumping of energy into residual modes. Success- 
ful application of this method of active damping 
requires that the control force locations be 
chosen so as to minimize residual energy. This 
point is discussed in more detail subsequently. 

To describe force apportioning mathematical- 
ly, we write the equations of motion for a linear 
structure discretized to n degrees of freedom 
(d.o.f.), 

( 1 ) 


Internal passive structural damping is omitted to 
simplify the analysis. ° We denote a target mode 
as mode s, and the vector of coherently phased 
actions (forces and moments) which selectively 
affects mode s as 0;(t). Substitiiting into Eq. (1) 
Qs and the standarS normal mode transformation 

n 

g(t) • [♦] C(t) Cp(t) (2) 

r-1 

leads, with mass and stiffness orthogonality con- 
ditions, to the modal equations of motion, 

• (♦piOj). r • 1,2,...,n (3) 

The mathematical objective of force apportioning 
is to select the constant shape fc of Qc(t) so 
as to maximize generalized force’(|;,„9s) into the 
target mode while minimizing generalized forces 
(ir*Q$)> r / s, into all residual modes, thus 
minimizing residual responses. Er(t)', r /t s, 
Existing methods5-8 for calculating fs have been 
developed primarily for application with experi- 
mental transfer function data. For purposes of 
this study, however, we present in Section V a 
simpler method which requires a priori knowledge 
of modal parameters. 

In the next section, expressions for Qs{t) 
are developed (with fs left unspecified) leading 
to standard viscous and Coulomb types of modal 
damping. 


III. ARTIFICIAL VISCOUS AND COULOMB 
TYPES OF DAMPING 

From Eq. (2), the portion of response in 
any particular d.o.f. 1 due to mode s is 

‘’is^‘> ‘ ♦is ^s^‘> 

We specify that motion qi(t) of a single selected 
translation, the ith d.o.f., is measured by a 
sensor, typically an accelerometer, and that the 
quantity qis(t) for each target mode is produced 
from q<(t) by narrow-band filtering around the 
target mode natural frequencies. 


III. A Artificial Viscous Damping 

To produce simple viscous type damping in a 
single target mode, we specify 

9s “ fs ^is (5) 

where Ks is a feedback constant to be determined. 
Thus, modal equations (3) become 


•Note, however, that the actual structural damping 
of large space structures is likely to be slight. 
If so, it will produce negligible coupling of the 
undamped normal modes® and, therefore, can only 
complement, rather than detract from, the damping 
provided by artificial means. Thus, the neglect 
of natural damping in the analysis should, for 
large space structures, lead to conservative 
results. 


[in] q + [k] q = Q 


sgn(.) 


ej, (6) 


r « 1 ,2,. . . iH 

In tht numerical simulations of Section VI, 
non-zoro initial displacement and zero initial 
velocity are specified for each target mode, and 
it is required that each target mode exponential 
decay envelope be damped to a portion Ps « t of 
its initial value in time tj> Thus, the required 
viscous damping factor is 




0 ) 


and, from the sth equation of (6), the feedback 
constant is 


,, 2 "s "i 


( 8 ) 


If a number m of target modes are to be 
damped simultaneously, then the appropriate 
action vector is 


9(t) 


III 

■E 

s»l 


*^s ti “lis 


(9) 


In Eq. (9), we assume for simplicity that the 
m target modes are the first m modes of the 
structure, but in general the target modes need 
not be only the lowest modes, nor need they be a 
series of consecutive modes. With Eq. (9), the 
modal equations become 


fMjii ^ 

r r r r r 


m 


Ks(9,.*/s^ *is «s* 
r » l,2,..,n 


where Ks for each target mode is calculated from 
Eqs. (7) and (8), Figure 1* is a conceptual 
block diagram illustrating the implementation 
of Eq. (10) for active damping of two modes of 
a discretized beam, with the use of three control 
thrusters. For simultaneous damping of more 
than one mode, as in Fig. 1, the elements of 
apportioning vectors fs, s = l,2,...,m, can be 
represented as elements of a rectangular matrix. 


1 for {.) » Cy 

0 for |{.)1 < Cy 

[-1 for (.) < -c„ 


and Cy is a threshold velocity level, which is 
zero for pure Coulomb damping, Small cv ^ 0 is 
considered here because any actual propulsive 
device must shut off at some non-zero level of 
command signal. Equation (11) is an approxima- 
tion to the type of propulsive thrust often 
referred to as "bang-bang" (or "on-off") thrust 
with a deadband^, with Eqs. (4) and (11), modal 
equations (3) become 

, . Kj(4^,fj) sgn(*|j55),(12) 


r » 1,2,. ..,n 


With the initial value of target mode dis- 
placement in the ith d.o.f, specified as 
dis(O) • 5s(0)/f1s *od zero initial velocity, we 
require that mode s be suppressed in time ts. 
Provided that Cy is small relative to peak target 
mode velocities, the feedback constant K$ which 
will give suppression time t$ is calculated as 
follows. For pure Coulomb damping (cy • 0). 
decay per period 2ir/we of 5e(t) in the sth equa- 
tion of (12) islO 

4 ♦is 

Ms 


The initial amplitude Us(0)| equals the decay 
per period times the number of periods to sup- 
pression, usts/2it. Thus, a little algebra gives 


S’ tj 


(13) 


If a number m of target modes are to be 
damped simultaneously, then' the appropriate 
action vector is 

m 

9(t) • Kjfj sgniq^j) (14) 
s=l 


Artificial viscous damping requires that 
the control forces be continuously variable in 
time, and it leads to a mathematically linear 
problem. It is possible, though, that propulsive 
control units will not be capable of delivering 
continuously variable thrust. So we consider 
next a type of control force which perhaps is 
more likely to be achieved in practice, but 
leads to a nonlinear problem. 

III.B Artificial Coulomb Oainpinq 


and the modal equations become 
m 

Kj($^,f5) sgn(^^5g,(lS) 

s*l 

r » 1,2,. ..,n 

where Ks for each target mode is calculated from 
Eq. (13K 


To produce damping similar to Coulomb 
dampinglO in a single target mode, we specify 

t 

9s * ■ fs ®9"^9js) (11) 

where 


* Figures and tables are located at the end of the 
paper, 


It should be noted that for small cy the 
individual control forces for single-mode damping, 
Eq. (11), will vary with time in nearly a square 
wave. However, because a number of nearly square 
waves are being superimposed in Eq. (14), the 
time variation of each control force for multi - 
mode damping generally will have the character 
of erratic stairsteps. This is illustrated in 
Section VI, 
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IV. SPECIAL CASE: VISCOUS DAMPING 
OF ONE TARGET MODE 

This rtiativily slmpU case is amenable, in 
part, to c1osed»form solution. The solution ob« 
talned will provide useful Information which Is 
applicable, at least qualitatively, to the analy- 
sis of more complicated cases. 

We seek solutions of modal equations (6) 
with feedback constant (8). The sth modal equa- 
tion Is rewritten as 

and we specify Initial conditions ($(0) ■ ^so M 
and ts(0) ■ 0. The solution Is 

^(t) “ So^"s^“dS^® * C0S(Wyjt + 9j)(16) 

where 

“ds ’ “s V’ ■ 

dj - tan'^l-CjUj/w^j) 

Next we substitute Eq. (16) Into the right-hand 
sides of all residual (r f s) modal equations of 
(6) and, again using Eq. (6), find the non-homo- 
geneous equations 

•Cut 

where 

" ' Vi~V 

With Initial conditions 5^(0) ’ 0 "“nd 5r(0) ■ 0 
for r / s, the solution of Eq. (17) can be written 
as a convolution integral, 

’ FtF’ /" ® * s1nu^5TSlnu^(t-t)dT(18) 

r 

Equations (16) and (17) show that the target 
mode eventually decays and excitation of the 
residual modes eventually ceases. Sut after the 
target mode Is suppressed, the residual modes 
remain In a steady state of free vibration. This 
steady-state motion can be determined by setting 
the upper limit In Integral (18) to Infinity, 
which also makes' the Integral reasonably tractable. 
The result of this evaluation has the form 

5/®(t) • 5^ cos(w^t + 8^) (ig) 

where ?r ^nd Sp are algebraically complicated. 

Rather than examining amplitude ?p. It is 
perhaps more useful to evaluate the residual 
energy, which Is a significant global quantity. 

The Instantaneous energy In any mode is 

Ep(t) = j Mp(jp^ + ip^) (20) 


Thus, the ratio of steady-state residual mode 
energy to Initial energy Is found to be 



Equation (21) demonstrates that selection of an 
active damping factor Involves a compromise! 
a higher ;$ li* desirable to suppress the target 
mode more rap*id1y, but a lower ;$ Is desirable to 
minimize residua] energy. Equation (21) also 
shows that residual modes with frequencies some- 
what greater than that of the target mode are 
relatively Immune to excitation, the residual 
energy varying nearly as the fourth power of the 
frequency ratio. On the other hand, residual 
modes with frequencies lower than that of the 
target mode are not protected by a favorable fre- 
quency ratio and seem somewhat susceptible. In 
the limit as wp 9oes to ze»’o, Eq. (21) gives the 
residual energy remaining In a rigid body mode. 



Equation (22) shows that force apportioning vector 
fs must be orthogonal to all rigid body mode 
shapes In order not to disturb the structure's 
rigid body orientation. 


V. FORCE APPORTIONING BASED ON A 
CONDENSATION TECHNIQUE 

We first establish an apportioning of forces 
and moments which Is applicable t) the full dis- 
cretized structural model characterized by Eq. (1). 
To Isolate the sth mode from all others, we set 

9s “ W tj 9(t) 

where g(t) Is an arbitrary function of time. Then 
modal equations (3) become 

Mp5p *■ M^Wp^ 5p » Mj g(t), r » l,2,...n 

where <Spj Is the Kronecker delta. Thus, the Ideal 
apportioning of forces and moments Is 

fs » Cm] 4s (23) 

But Eq. (23) Implies the presence of a control 
force-or moment at each d.o.f, of the structure, 
which is generally Impossible. So we develop 
next a condensed, or reduced order, mathematical 
model of the structure, for which an apportioning 
analogous to Eq. (23) can be defined. 

The condensed model Is required to include 
nc < n translation d.o.f. corresponding to control 


4 


thrustirs ind to preserve wUhout alteration nc 
selected /nodes of the structure. Accordingly, 
we define the nc x nc condensed inertia and 
stiffness matrices by analogy with the standard 
orthogonality conditions, 

where [•>‘^‘'1 Is the nc x nc partition of the full 
modal matrix [♦] containing only the d.o.f. and 
modes to be retained In the condensed model, and 
C*MC J and C“u 2 mc j are the nr x nc partitions of 
the generalized mass and stiffness matrices, re- 
spectively containing only parameters of the re- 
tained modes. Provided that not singular, 

the condensed Inertia matrix is therefore 

C")“} * (25) 

and Is calculated similarly. 

Thus, the condensed version of Eq. (1) Is 

9“= - f (26) 

where g*^and are the partitions of g and Q, re- 
spectively, which contain only the translation 
d.o.f, and the associated forces retained In the 
condensed model. It can be provedH that problem 
(26) preserves exactly the modal parsnieters and 
the free and forced response of the retained modes 
and d.o.f. 

In analogy with Eq. (23), we define a con- 
densed force apportioning vector which completely 
isolates mode s from all other modes retained In 
the condensed model , 

fs“ * [m'==] is" (27) 

To simplify Eq. (27), we define in Eq. (25) 

Cv] * or (j^, (28) 

Substituting Eqs. (25) and (28) Into (27) then 
gives 

fs“ = «S ^s <29) 

The constant Mj in Eq, (29) Is Irrelevant since. 

In calculations, fc<^c i$. normalized to a maximum 
element value of tl, as are all mode shape vectors. 
The nc X 1 vector jj the non-zero partition 
of the full n X 1 apportioning vector fs, all 
remaining elements of fj being set to zero. Since 
(ir>fs) * Msips for all modes retained In the 
condensed model,' this apportioning will control 
the target mode without affecting the other re- 
tained modes. But, In general, It will excite 
the modes which have not been Included In the con- 
densed model. When used in Eqs. (10) or (15), 
this apportioning completely controls ny selected 
modes with Oq selected control thrusters (pro- 
vided that Is not singular), but all other 

residual modes are uncontrolled. An Important 
objective in control thruster selection, therefore. 
Is to minimize excitation of uncontrolled modes. 


The condensation-apportioning procedure de- 
scribed In this section requires complete Informa- 
tion on the modes to be retained In the reduced 
order model, and It assures that each of the re- 
tained iiiMdes can be Isolated from all other re- 
tainec! .roles. On the other hand, apportioning 
ttchniO',e- used In modal testing®'® do not require 
nioi.al i.’nmeters a priori . but they also do not 
pro iur.e /ny assured measure oV mode Isolation. So 
Eq. (29) probably represents the mathematically 
best apportioning of a limited number of control 
thrusters, but It may not be best In practice be- 
cause It requires previously determined modal 
parameters. 


VI. NUMERICAL SIMULATIONS OF ACTIVE DAMPING 

Time history simulations have been calculated 
for evaluation of active damping by force appor- 
tioning, the performance measure of interest being 
the residual energy remaining after the target 
modes are suppressed. 

Figure 2 depicts the unrestrained plane gril- 
lage structural model used in the simulations, It 
is a relatively simple mathematical model de- 
signed to exhibit structural dynamic character- 
istics typical of proposed solar power sate 111 tesl 2 ^ 
in particular, low natural frequencies and high 
modal density. The model Is an assemblage of 
straight bending-torsion beam finite elements with 
twelve nodes, as numbered In Fig. 2, Each node 
has a translational d.o.f. In the Z direction and 
rotational d.o.f, about the X and Y directions, 
for a total of 36 d.o.f. Accordingly, the model 
has a rigid body translation mode In the Z direc- 
tion and rigid body rotation modes about the X 
and Y directions. Each of the six square bays 
has side length L = 4000 m. A mass of 4,8 x 10® 
kg (lot of the total mass of beams) with rotational 
Inertias of 3.0 x lOlO kg • m^ is concentrated 
at node 2. Inertias of the beam elements are 
represented by consistent mass matrices. Each 
element has mass per unit length of 35.3 kg/m, tor- 
sional inertia per unit length of 5,65 x lo4 
kg > m, bending stiffness El of 7,55 x 10l2 N • m^, 
and torsional stiffness GJ of 5.81 x I 0 l 2 N ■ m2. 
Natural ,‘requencies and periods of the first eleven 
flexible modes are listed in Table 1. Mode shapes 
of the first six flexible modes are shown In Figs. 3. 

For each simulation, the target modes selected 
were modes 4-6, the first three flexible modes. 
The only non-zero Initial conditions specified 
were mode 4 Initial value 54 ( 0 ) corresponding to 
node 1 translation qiqlO) * 1.0 m, initial 
values of 55 ( 0 ) and 56 ( 0 ) corresponding to equal 
initial potential energy in each target mode. The 
feedback signal chosen was node 1 translation 
velocity qi(t). In addition to the target modes 
4-6, rigid body modes 1 - 3 were retained in 
the condensed model In order to avoid disturbing 
the rigid body orientaion (cf. Eq. (22)). Six 
thrusters were required to control these six modes. 

In calculations, only the first 14 of the 
total 36 modes were retained, i.e. Eq. (2) was 
truncated after r » 14. .This truncation is 
justified by the (us/up)‘l term in Eq. (21) and by 
the observation that the highest retained modes 
contributed negligibly to the total response. The 
modal equations of motion were integrated by a 


5 


fourth-order Runge-Kutta method with constant step 
size (2 s for viscous damping, 1 s for Coulomb). 

Active damping performance was evaluated for 
different control thruster locations, different 
times to damp the target modes, and different 
damping types, viscous or Coulomb. Case 1 Is the 
reference case, and the parameters were varied 
Individually for Cases 2-4. In Case 1, control 
thrusters are located at nodes 1, 2, 3, S, 10, 
and 12, and damping Is viscous with time t$ > 

400 s to damp each target mode down to portion 
Ps « 0.01 of Its Initial value. Time histories 
for selected quantities are shown in Figs. 4. 
Individual modal energy values were calculated 
from Eq. (20) and summed to produce Fig. 4a, the 
history of total energy. The active damping In 
this case Is very effective, with residual energy 
being less than 3% of Initial energy. 

For the damping type and time of Case 1, 
several combinations of six control thruster loca- 
tions were Investigated. A few combinations can- 
not be used because tlie associated C^cc] matrices 
are singular (see Eqs. (28) and (29)), but most 
can be usedli. The thruster locations of Case 1 
produced one of the best results. For Case 2, 
a control thruste^r was moved from node 5 to node 
11, and this change produced one of the worst 
results. Figures 5 show that the active damping 
Is counterproductive; the residual energy Is 28 
times the Initial energy, with. b 73S of the residual 
energy remaining In mode 8 alone. This poor per- 
formance can be explained by examination of the 
force apportioning vectors and feedback constants 
for Cases 1 and 2 (Table 2) and the mode shapes 
for modes 4 and 8 (Figs. 3 a,e). Equations (10) 
and (21) suggest that ('('s.fs) should be large in 
comparison with (■^rtfs)', r s, For target mode 
4 and residual mo'de 8 in Case 1, this condition 
Is satisfied. In Case 2, however, (M.fq) Is 
very small In comparison with (48»-f4) ®nd, con- 
sequently, K 4 Is relatively large and a great 
deal of energy is pumped into mode 8. 

The parameters of Case 3 are those of Case 1, 
except that damping time Is halved, t«. * 200 s, 
so that Cs is doubled. Although Eq. (21) sug- 
gests that the residual energy might be quadrupled 
by this Increased damping, the actual residual 
energy in Fig. 6a is only a little more than 
twice that of Case 1 . 

The parameters of Case 4 are those of Case 1 
except that damping Is the Coulomb type, Eq. (14), 
with --V = 0,0001 m/s for all target modes. 

Figures 7 demonstrate that this damping Is very 
effective. Comparison of Figs. 4c and 7c shows 
that lower force levels are required for Coulomb 
type damping, b.ut that these levels must be 
maintained throughout the damping period. For 
t •> 400 s. Fig, 7c exhibits thruster chattering, 
which results from one target mode alternating 
between free and damped motion around the thresh- 
old velocity level . 


VII. CONCLUDING REMARKS 

The basic form of active vibration damping 
by force apportioning has been described with 
standard structural modal analysis and numerical 


simulations. For possible application to satel- 
lite structures, rigid body modes have been con- 
sidered, but the Influence of gravity on orbital 
position and attitude has been neglected. Feed- 
back control laws producing modal damping of the 
simple viscous and Coulomb types have' been used. 

To permit evaluation of the method In Its 
basic form, only Idealized conditions have been 
considered, It has been assumed that all neces- 
sary modal parameters are known and that sensors, 
actuators, and narrow-band filters have Ideal 
operating characteristics. The damping of two or 
rore nades with effectively Identical natural fre- 
quencies has not been Investigated. 

Under these Idealized conditions, the method 
Is effective and relatively simple. It uses a 
single sensor and a limited number of discrete 
thrusters to control completely an equal number 
of modes. The significant measure of performance 
Is the amount of residual energy fed Into the un- 
controlled modes during the process of damping 
the controlled modes. In the numerical example 
presented, a judicious choice of control thruster 
locations and moderate damping of the controlled 
modes produces only an Insignificant amount of 
residual energy. 
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Table 1 Flexible mode natural frequencies and 
periods for unrestrained plane grillage 
structure 


Mode 


Frequency (Hz 


4 

0,0062 

161,29 

5 

0.0065 

153.85 

6 

0,0129 

77,52 

7 

0.0143 

69.93 

8 

0.0155 

64.52 

9 

0.0180 

55,56 

10 

0,0217 

46.08 

11 

0.0232 

43.10 

12 

0.0245 

40.82 

13 

0.0448 

22.32 

14 

0.0489 

20.45 

Table 2 

Force apportioning vectors and feedba 
constants for active damping of mode 


Case 1 

Case 2 

^1.4 

-0.128 

-0.244 

^2.4 

0.923 

0.488 

^3.4 

-0.128 

-0.244 

'5,4 

-1.000 

- 

'l0,4 

0.167 

0.500 

'll. 4 

- 

-1.000 

^2.4 

0.167 

0.500 

K.(kg/sec) 5.40 x 10^ 

1.54 X 10® 
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Figs. 4 Case 1: control thrusters 1, 2, 3, 5, 
10, 12; viscous damping with tc * 400 s, 
Ps * 0.01 


Figs. 5 Case 2: control thrusters 1, 2, 3, 10 
11, 12; viscous damping with te * 400 
Ps • 0.01 
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APPENDIX B 


B.l IDEAL FORCE APPORTIONING 

(This section Is Intended to be the replacement for Section V of Ref, 2. 
Accordingly* the section title number* equation numbers* and reference num- 
bers are those of Ref. 2.) 


V. IDEAL FORCE APPORTIONING 

We define “controlled" modes as being the m modes described by the Index 
s In Eqs. (10) and (15). Although this set has been defined previously as 
consisting only of the target modes* It Is clear now* In view of Eq, (22)* 
that the set of m modes should include not only the flexible target modes* 
but also all appropriate rigid body modes. It is assumed in the following 
that the m controlled modes are specified. 

In order to establish a mathematically Ideal force apportioning vector 
fs for a single controlled mode (denoted by subscript s)* we define Qs(t) » 

M Ml 

f$ g(Cs)i where g(c$) represents the feedback time dependence descr'Ibed In 

Eqs. (5) and (11). It Is generally possible to provide control forces at 

only a few translation d.o.f. of a discretized structure. Hence* the n x 1 

vector fg must usually Include mostly zero elements. We therefore denote 

the actual number of control thrusters (actuators) positioned on the structure 

a 

as ng, and the corresponding n x 1 force apportioning vector as f * the 

** ~s 

elements of which should generally be non-zero. Accordingly, the n^ x 1 
partition of modal vector corresponding to the control actuator d.o.f. Is 
denoted Therefore* Eq. (3) becomes 

= (*J* fg) g(55), r = 1, 2, n (23) 

It Is necessary in the present context to choose the number of control 
actuators equal to the number of modes to be controlled, n, - m. Then the m 

a 


equations of (23) which describe the controlled inodes can be written In 
matrix form as 

g(L) (24) 

W t 9 

Where and are appropriate m x m partitions of the n x n 
generalized mass, generalized stiffness, and modal matrices, respectively. 

To decouple controlled mode s from all other controlled modes In Eq. (24), 
we seek f^ such that only the modal equation for mode s has a non-zero right- 
hand side. Clearly, this objective Is achieved by setting 

.1 ■ <«) 

where ^^1s the column vector corresponding to control led modes of the matrix 

Ct] • (26) 

In calculations, f^ Is normalized to a maximum element absolute value of 1, 
as are all mode shape vectors. Provided that Is not singular, then the 
entire apportioning matrix of Eq. (26) can be calculated; then the set of 
vectors f^, s ° 1, 2, ... m, from Eq, (25) decouples the m equations of the 
controlled modes In modal Eqs. (10) or (15), but leaves the remaining n-m 
equations of the uncontrolled (residual) modes excited by the controlled 
modes. An Important objective In control thruster selection, therefore. Is 
to minimize excitation of the uncontrolled modes. 

The apportioning procedure described In this section requires complete 
modal Information on the controlled modes, and It assures that each of the 
controlled modes can be isolated from all other controlled modes. On the 

5 8 

other hand, apportioning techniques used in modal testing do not require 
modal parameters a priori , but they also do not produce any assured measure 
of mode isolation. So the apportioning calculated by Eqs. (25) and (26) 


23 
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represents a nia thematically Ideal apportioning of a limited number of control 
thrusters, but It may not be best In practice because It requires previously 
determined modal parameters. 

B.2 NON- IDEAL FORCE APPORTIONING BY MATRIX PSEUDO- INVERSE 

We begin with a somewhat more explicit derivation of the results of Eqs. 
(25) and (26) above. The complete m x m force apportioning matrix Is defined 
to be 

tF’] - [ff fj] (B.2.1) 

To uncouple ell controlled Diodes, we seek [f‘] such that (see Eq. (24) above) 
[F“] " Cy (B.2. 2) 

Where Is the m x m identity matrix. Clearly, 

[F“] • (B.2. 3) 

In our computations for fj, s « 1, 2, ...» m, each column of [F®] Is subse- 
quently normalized to a maximum element absolute value of 1. 

In practice, It will probably be desirable to control a greater number 
of modes than there are control actuators, m > n.,. For this situation, we 

a 

define the force apportioning matrix [F^] just as In Eq. (B.2.1), except that 
now each column has length n^. We again seek the solution of Eq. (B.2. 2), 
where now Is an n^ x m matrix. Since the unknowns in this equation are 
overdetennined, a unique, exact solution does not generally exist. We pro- 
ceed formally as follows: 
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(B.2.4) 


c*“] [♦“f Cf*] ■ [♦“=] Dy 

CF*1 - ([♦*'] C4*"=]^)'' [♦“] 

Error jj Is defined es the jth column of [F°] - Then Eq. 

(B.2.4) represents the least-squares approximate solution of Eq, (B.2.2) for 
m > n^j, for which each norm (e^, gj), j » 1, 2, ..., m, Is a minimum.'*' 

Eq, (B.2.4) clearly reduces to Eq. (B.2.3) If m » n,. Eq. (B.2.4) Is a form 

a 

of the so-called matrix pseudo- Inverse for the overdetermined least-squares 
case, If It Is desired to weight the control of certain modes over control 
of others, a generalization of Eq. (B.2.4) for a weighted least-squares ap- 

it 

proximate solution is easily developed. 


•k' 

W. L. Brogan, Modern Control Theory , pp. 90-92, Quantum Publishers, 1974. 
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APPENDIX C 

ACTIVE DAMPING OF HIGHER-FREQUENCY MODES 

A single case was Investigated, the numerical simulation model being 
the unrestralnev plane grill age described In Ref. 2 . The target modes were 
specified to be modes 7, 8, and 9 (rather than modes 4, 5, and 6. as In Ref. 2), 
and the controlled modes were specified to be the target modes plus the rigid 
body modes, modes 1, 2, and 3. 

Twenty-five combinations of six control thruster locations were consi- 
dered in a trial -and-error search for the combination producing the least 
excitation of the uncontrolled modes. The combination of thrusters I, 2, 3, 

5, 10; 12, which Is quite effective If nodes 4, 5, and 6 are the target modes 
(cf. Ref. 2), was found to be very poor for target modes 7, 8, and 9. The 
most effective combination for this case (among those Investigated) consists 
of control thrusters 1, 2, 3, 7, 9, 11 (cf. Fig. 2 of Ref. 2), and this 
combination was used In time history numerical simulations of active damping. 

Viscous-type damping was used. The Initial conditions specified cor- 
respond to an equal amount of initial potential energy In each of the target 
modes, with zero initial kinetic energy In the target modes, and zero initial 
energy In all non-target modes. The feedback signal used was node 1 transla- 
tion velocity q^(t). The exact force apportioning vectors (Eq. (29) of Ref. 2) 
were used, and filtering was taken to be Ideal. 

Figures C.l and C.2 are energy-decay time histories comparable to Figs. 
4a-7a of Ref. 2. With specified time t^ = 400 s to damp each target mode 
down to portion p^ = 0.01 of its Initial amplitude, the residual energy shown 
on Fig. C.l is about 21% of the initial energy. About half of the residual 
energy resides in each of modes 6 and 10, the modes immediately adjacent to 
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the target n»des. Figure C.2 shows that the doubling of specified damping 
time to 800 s produces a much more satisfactory final residual energy of 6 % 
of the Initial energy. 
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APPENDIX D 

INEXACT NODAL PARAMETERS FOR THE CONTROLLED STRUCTURE 

We have Investigated this subject pritnarlly by doing numerical experiments » 
I.e.t computer simulations of. active damping* rather than a theoretical study. 

Our first step was to define a ''model'' structure and an "actual" structure, 
both of which are really just mathematical models. The "errors" In the "model" 
structure relative to the "actual" structure are small differences In the mass 
distributions. These "errors" are presumably typical of the Inaccuracies 
which Inevitably exist In a finite element model relative to the real structure. 

Following definition of the "model" and "actual" structures and calcula* 
tion of the. vibration modes for both, the Ideal force apportioning vectors 
were calculated by Eqs. (28) and (29) of Ref. 2 from the modal matrix of the 
" model " structure. Finally, these force apportioning vectors were used In 
numerical simulations of active damping on the " actual " structure. The quality 
of the active damping then was an Indication of the effect of non-ideal force 
apportioning, which results from Inexact modal parameters. A short digression 
on the theory is appropriate before presentation of numerical results. 

The mathematical consequence of non-ideal force apportioning Is to pro- 
duce non-zero Inner products (i^, f^) for r s In the modal equations of the 
controlled modes (those modes denoted by the Index s In Eqs. (10) and (15) of 
Ref. 2), thus coupling the controlled modes. It Is to be expected, then, that 
active damping with non-ideal force apportioning will be less effective than 
active damping with Ideal force apportioning. Moreover, the coupling of con- 
trolled modes could conceivably lead to their Instability, which Is not pos- 
sible without coupling. Nevertheless, most of the numerical simulations to 
be described next exhibit neither substantial damping degradation nor Instability. 
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The "actual" structure considered here (referred to as Model 5} Is 
slightly different than the structure described In Ref. 2 (referred to as 
Model 3)» for reasons to be given below. Model 5 Is Identical to Model 3 In 
geometry and mass distrlbutloni and Its stiffness distribution Is the same In 
all respects except that Its center longitudinal member (conne Mng nodes 
2* 5, 8, and 11 In Fig. 2 of Ref. 2) has a greater torsional stiffness, 

6J ■ 2.0 X 10^^ N.m^. This stiffening separates mode 5, the first flexible 
torsion mode, from mode 4, the first flexible bending mode (cf. Figs. 3 and 
Table 1 of Ref. 2); mode 5 natural frequency of Model 5 Is 0.0097 Kz, In 
contrast to 0.0065 Hz for Model 3. The frequencies of bending modes 4 and 
6 are Identical for the two models. 

The "model" structure considered Is Identical to the "actual" structure 

except that It has additional small lumped masses at nodes 8, 10, and 11* 

(cf. Fig. 2 of Ref. 2)i each has mass of 2.0 x 10^ kg and rotational Inertias 
9 2 

of 8.0 X 10 kg»m . The effects of these mass perturbations are to make the 
"model" structure's lowest natural frequencies on the order of 10% lower than 
those of the "actual" structure, and also to alter the mode shapes somewhat. 

Ideal force apportioning vectors were calculated for both the "actual" 
and "model" structures, with the controlled modes specified as modes 1-6 and 
the control thruster locations specified as nodes 1, 2, 3, 5, 10, and 12. 

These vectors are listed for target modes 4, 5, and 6 In Table 0.1. Com- 
parison of the vectors for "actual" and "model" structures shows that those 
for bendlhg modes 4 and 6 are quite similar, while those for torsion mode 5 
differ substantially. 

* 

It would be more natural tf the "model" structure were symmetric and the 
"actual" structure had the asymmetric mass perturbations, but It proved 
convenient In computations to proceed as described above. 


30 


Energy time histories for viscous and Coulomb active damping are shown In 
Figs. D.1 and 0,2, respectively. Simulations for both Ideal and non-ideal 
apportlonings were run to provide a comparison. Viscous damping with the non- 
ideal apportioning Is, surprisingly, slightly better than that with Ideal 
apportioning. Coulomb damping with non-ideal apportioning produces a slightly' 
faster decay but also excites slightly more residual energy than that with 
Ideal apportioning. 

The results In Figs. D.1 and D.2 are representative of a number of similar 
cases Investigated In which the modes of the "model” structure are not sub- 
stantially different from those of the "actual" structure. The Invariable 
observation was that non-ideal apportioning resulting from Inexact modal 
parameters produces a slightly different damping character than does Ideal 
apportioning, but never any serious degradation of short-term damping quality 
nor any noticeable long-term Instability. The conclusion then is that the 
quality of active damping by force apportioning Is not very sensitive to small 
errors In modal parameters of the controlled modes. 

However, It Is not valid to state that active damping quality Is always 
Insensitive to small errors In the primitive physical parameters (I.e. , stiff- 
ness and mass distributions) of the controlled structure; the existence of 
closely spaced modes (high modal density) can cause problems. To illustrate 
this, we evaluated a case for which Model 3 Is the "actual" structure, and the 
"model" structure differs only In the presence of a small lumped mass at node 10. 
An Important effect of this small additional mass Is to reverse the characters 
of the two closely spaced, lowest flexible modes of the "model" structure 
relative to the same modes of the "actual" structure (cf. Table 1 and Figs. 3 
of Ref. 2): mode 4 of the "model" structure Is torsion, and mode 5 Is bending. 
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Because like-numbered mode shapes of the "actual" and "model" structures are 
fundamentally different In this case, active damping of the "actual" structure 
by force apportioning calculated for the "model" structure produces very poor 
results. For this particular example, the problem can be solved simply by 
reversing the order of mode shapes 4 and 5 In the modal matrix of the "model" 
structure before calculating the force apportioning vectors.* In general, 
however, the existence of substantial errors In estimates of mode shapes can 
seriously degrade the quality of active damping by force apportioning. It 
seems probable that actual hardware Implementation of this or any other scheme 
for active vibration control of large space structures will require an adaptive 
capability. 


* 

Nevertheless, we considered It desirable to avoid this confusing complication 
by using Model 5, which has well separated first bending and torsion modes, 
for the principal numerical examples of this appendix. 


Table D. 1 Ideal Force Apportionings 
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"Actual" 
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Structure Hode 
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"Model" 
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Structure Node 
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0.15 
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0.59 
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-1 

0.95 
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-1 
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-0.13 

-1 

0.57 

-0.11 

-1 

0.54 

5 

-1 
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-0.21 

-1 

-0.07 

-0.18 

10 

0.17 

-1 

0.03 

0.20 

-0.56 

0.01 
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(a) Force apportioning for the "actual" structure (b) Force apportioning for the "model" structure 


Fig. D.l 


The effect of non-ideal force apportioning on damping of the "actual" structure: Model 5 with 
viscous damping; controlled inodes 1-6. target modes 4-6; control thrusters 1, 2. 3, 5, 10, 12; 





0.0 




(a) Force apportioning for the "actual" structure (b) Force apportioning for the "model" structure 


Fig. D.2 The effect of non-ideal force apportioning on damping of the "actual" structure: Model 5 with 
Coulomb damping; controlled modes 1-6, target modes 4-6; control thrusters 1, 2, 3, 5, 10, 12; 
equal initial potential energy in modes 4, 5, and 6; t = 4C0 sec, = 0.0001 m/sec 
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APPENDIX E 

NON-IDEAL' FREQUENCY FILTERING 

Equations (ID) and (11) of Ref. 2 can bt written as 

+ E S<*r* 

s ■ 1 


1i 2| • • • t ^ 


(E.1) 


where 


h(qis) 


for viscous damping 
sgn(q^$) for Coulomb damping 


The term q^^ In Eq. (E.1) Implies Ideal filtering at each controlled mode 
natural frequency s ■ 1» 2, ...» m» as Is Illustrated conceptually on 
Fig. 1 of Ref. 2. The Ideal filter for controlled mode s passes only the 
q^^ spectral component of the total sensor velocity, q^ > The 

Ideal filter, then, has an Infinitely narrow pass-bank and Infinitely steep 
rolloff rates. Since q^^ ■ a mathematical characteristic of Ideal 

filtering Is that It does not couple the controlled modes In Eq. (E.1). 

Those modes might be coupled If non-ideal force apportioning Is used with 
Ideal filtering, but In such a case the non-ideal apportioning Is responsible 
for the coupling. If both apportioning and filtering are Ideal, then the 
controlled modes are not coupled at all. 

On the other hand, a non-ideal, more realistic filter has a finite 
pass-band and finite rolloff rates, and It consequently passes spectral com- 
ponents- of all frequencies. The characteristics of a non-ideal filter for 
controlled mode s are Illustrated on Fig. E.1. For such a filter, the 
function h(q^^) In Eq. (E.1) must be replaced by a filter output function 

F,jg ” ^ ^ **'* 
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where constants A^j are attenuation factors Illustrated In Fig. £.1. For 
Ideal filtering* the . Kronecker delta. Calculation of the at- 

tenuation factors for non-ideal filtering Is discussed below. For non- 
Idoal filtering* then* the modal equations have the form 






r « 1* 2* ...* n (£.2) 

'n £q. (£.2)* both controlled and uncontrolled modes are coupled by the 
non-ideal filtering* regardless of the type of force apportioning used. 

So It is to be expected that non-ideal filtering will reduce the effective- 
ness of active damping by force apportioning and* moreover* that It may even 
produce system Instability due to the coupling. Nevertheless* numerical 
results presented below for typical non-ideal filter characteristics ex- 
hibit neither substantive degradation of damp1n<j quality nor instability. 

To calculate constants A^j* one must specify frequency cutoffs and 
rolloff rates. Filter rolloff rates are designated In terms of decibels 
per octave. So* for a high frequency rolloff at rate R db/oct* we have 
(see Fig. £.1) 

20 ■ - R 

» - R 1og^Q(uj/u5,^)/log^Q2 

Solving for the attenuation factor gives* 


- R/20 log..2 

*sj * ’ “4 “sh 

Similarly, for a low frequency rolloff*. 

''sj ' ‘‘’j - Sjl 


(E.3a) 


(E.3b) 
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Wt hivt run several numerical simulations of active damping by 
force apportioning with non-ideal filtering. The structural models used 
In these simulations are those described In Appendl.n D. Simulations with 
Model 3 tested the effect of non-ideal filtering alone* since Ideal force 
apportioning was used, and simulations with Model S (as the "actual” 
structure) and the related "model" structure tested the effects of non- 
ideal filtering and oon-ldeal apportioning acting simultaneously. In all 
cases, modes 1-6 were specified as the controlled modes and modes 4-6 as 
the target modes, and control actuator locations were specified as nodes 
1, 2, 3, S, 10, and 12. 

The filter cutoff frequencies specified for both Model 3 and Model 5 
can be described with reference to Fig. E.2. A low-pass filter was used 
for mode 4, and band-pass filters were used for modes 5 and 6. The cut- 
off frequencies were selected to be 

“4h “ “5 j2 " ^“4 * “5^/^ 

“5h * “il ■ <“6 * 

“6h ■ ^“6 * “7^/^ 

The particularly Interesting characteristic of Model 3 Is the close proximity 
of natural frequencies 4 and 5 (cf. Table 1 of Ref. 2). This closeness 
might lead one to expect that non-ideal filtering would produce substantial 
coupling between modes 4 and 5 and, hence, would reduce the quality of 
damping. 

To test this expectation, we ran Model 3 cases with filter rolloff 
rates varying from the unrealistically high value of 160 db/oct, through 
the practically realistic values of 80 and 48 db/oct, and down to the un- 
realistically low value of 10 db/oct. Cases were run with both viscous 
and Coulomb types of damping, and with different Initial conditions. 
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Tht adverse effect of non-ideal filtering was found to be relatively 
small In all Hodel 3 cases w1th> typical rolloff rates of 80 or 48 db/oct. 

The energy time history for the very ivorst of these cases Investigated Is 
shown In Fig. E.3, which may be compared with Fig. 4a of Ref. 2. Whereas 
Ideal filtering was used for the latter, non-ideal filtering with 48 db/oct 
rolloff rate was used for the former; however, all other conditions were 
Identical for the two runs. Comparison of the two figures shows that non- 
ideal filtering In this case retards the rate of active damping, but does 
not Increase the final quantity of residual energy remaining In the structure. 

Another Model 3 case of Interest Is one for which the unrealistically 
low rolloff rate of 10 db/oct was used. The energy time history Is shown 
In Fig. E.4, which may be compared with Fig. 7a of Ref. 2. Ideal filtering 
was used for the latter and non-ideal filtering for the forfner, but all 
other conditions were Identical for the two runs. Figure E.4 shows a short- 
term rate of damping very similar to that of Fig. 7a of Ref. 2. But the 
residual energy Is greater In the case of non-ideal filtering; and. In 
fact, the energy Increases after reaching a short-term minimum (at about 
700 sec), strongly suggesting the presence of a system Instability. This 
Instability seems relatively weak and manifests Itself only after the In- 
itially existing vibrations are essentially suppressed, so It might be 
regarded as a “long-term" Instability. 

Numerical simulations of active damping with Model 5 tested the sit- 
uation* In which both the force apportioning and the filter cutoff frequencies* 
designed for the “model" structure were used for the “actual" structure. 

The filter cutoff frequencies were designed for the “model" structure In 
accordance with Eqs. (E.4), with one exception: the high frequency cutoff 
for the Mode 6 filter had to be adjusted upward slightly for use with the 
"actual" structure, because the computer program requires the natural fre- 
quency to be within thr^ pass-band. 
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Comparison of Figs. E,5a and E.Sb shows that In this case the presence of 
both non-ideal apportioning and non-ldeal filtering with a rolloff rate of 
80 db/oct only slightly reduces the quality of active damping. The case 
was repeated for rolloff rates of 48 and 24 db/oct, with almost no further 
deterioration of damping quality. One additional halving of the rolloff 
rate, to the unrealistically low value of 12 db/oct, resulted In the energy 
tima history of Fig. E.5c, which exhibits both substantative loss of short- 
term damping quality and long-term Instability. Finally, a Model 5 case 
was run for which the rolloff rate was set to 0 db/oct, i.e., no filtering 
at all. For this case, the entire system, consisting of both controlled 
and uncontrolled modes. Is strongly unstable. 

To sumtnarize and reiterate the observations made relative to non- 
ideal filtering In the particular cases studied: (1) the quality of active 

damping Is not substantively decreased by typical filter characteristics; 
(2) both poor damping quality and long-term instability may occur >r 
filtering with unusually low rolloff rates. 
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Fig. E.3 The effect of non-ideal filtering: Model 3 with viscous damping; non-ideal filters 

with 4S db/oct rolloff rate, cutoff frequencies (Hz) of 0.00638 for mode 4 low-pass, 
0.00638 and 0.00974 for mode 5 band-pass, 0.00974 and 0.01363 for mode 6 band-pass; 
controlled modes 1-6, target modes 4-6; control thrusters 1,2, 3, 5, 10, 12; equal 
initial potential energy in modes 4, 5, and 6; t^ = 400 sec, p^ = 0.01 
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Fig. E.4 The effect of non-ideal filtering: Model 3 with Coulomb damping; non-ideal filters 

with 10 db/oct rolloff rate, cutoff frequencies (Hz) of 0.00638 for mode 4 low-pass, 
0.00638 and 0.00974 for mode 5 band-pass, 0.00974 and 0.01363 for mode 6 band-pass; 
controlled modes 1-6, target modes 4-6; control thrusters 1, 2, 3, 5, 10, 12; equal 
initial potential energy in modes 4, 5, and 6; t = 400 sec, e = 0.0001 m/sec 
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(a) Ideal filtering, and force apportioning for the (b) Non-ideal filtering with 80 db/oct rolloff rate, 

"actual" structure and force apportioning for the "model" structure 


Fig. E.5 The effects of both non-ideal apportioning and non-ideal filtering on damping of the "actual" structure: 
Model 5 with viscous damping; filter cutoff frequencies (Hz) of 0.00713 for mode 4 low-pass, 0.00713 and 
0.01029 for mode 5 band-pass, 0.01029 and 0.01294 for mode 6 band-pass; controlled modes 1-6, target 
modes 4-6; control thrusters 1,2,3,5,10,12; equal initial potential energy in modes 4, 5, and 6; 
tg = 400 sec, p^ = 0.01 
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